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Abstract. We discuss the relation between multi-particle production in the Glasma 
at next-to-leading order and the physics of plasma instabilities. 



1. Introduction 

This talk is based on recent work on two topics: i) a formalism to compute multi-particle 
production in heavy ion collisions at next-to-leading order with Frangois Gelis [1] and 
ii) 3-l-l-D numerical simulations of Yang-Mills equations with Paul Romatschke [2] that 
investigate the role of plasma instabilities in the Glasma. The word "Glasma" describes 
the physics of the strongly interacting matter from the time when particles are produced 
in the shattering of two Color Glass Condensates (CGCs) in a heavy ion collision to 
the time when the Quark Gluon Plasma (QGP) thermalizes [3]. The theme here is that 
there exists a deep connection between topics i) and ii). It may provide insight into the 
puzzle of how the QGP is formed from the CGC at apparently very early times of < 1 
fm. 

2. Ab initio treatment of multi-particle production in heavy ion collisions 

Several thousand particles are produced in a heavy ion collision at RHIC and will soon 
be at the LHC. Can we understand how these particles are produced and follow their 
evolution and formation of a thermalized QGP from first principles ? Over a decade 
ago, heavy ion collisions were widely viewed to be an intrinsically non-perturbative 
problem not amenable to systematic computation of the properties of hot and dense 
matter formed. However, simple considerations have long suggested that the dynamics 
of partons at high energies are controlled by a semi-hard scale, the saturation scale 
Qs{x), which grows as x decreases thereby making a weak coupling description feasible: 
dsiQs) ^ 1- The physics describing the dynamics of partons at small x can be 
formulated as an effective field theory-the Color Glass Condensate (CGC) [1]- which has 
been applied to study particle production across the board among high energy scattering 
experiments. The basic ingredients in the CGC description of the nuclear wavefunction 
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at high energies consists of large x static color sources p(x^) and small x dynamical gauge 
fields. The sources are strong sources of strength p ~ l/^/os ^ 1. The probability of 
multi-particle production at high energies can be computed in the framework of particle 
production in field theories with strong external sources. Specifically, we find [1] 

p„=exp(-i5i:*„)ti Y. m 

where hr denotes the sum of vacuum-to-vacuum graphs with r cuts. This formula has 
remarkable features which are particular to field theories with strong sources: a) P„ is 
non-perturbative in g even for g <^ 1-no simple power expansion in terms of g exists, b) 
P„, for any n, gets contributions from cut tree vacuum graphs-this would not apply for 
field theories in the vacuum, c) Even at tree level, Pn is not a Poisson distribution. This 
last result is non-trivial because one might naively assume that classical field theories 
would have only "trivial" Poissonian correlations. The simple formula in Eq. [T]is all one 
needs to understand many features of the well known AGK calculus of multi-particle 
production [1]. 

Interestingly, even though the probabilities in Eq. [T]are non-perturbative for (7 -C 1, 
a systematic expansion in powers of g exists for moments of the multiplicity. The leading 
term for the first moment, the average multiplicity (n) = Ylm'^Pni is of order QiXjg^^ 
in the coupling but contains all orders in gp. The next-to-leading order contribution 
is ©(^f^) in the coupling but again all orders in gp. Another interesting feature of 
the moments is that even though they contain sums over r-particle Feynman graphs 
(corresponding to the in Eq. [1]), they can be represented simply in terms of solutions 
of equations of motion with retarded boundary conditions. At leading order, these are 
solutions to the Yang-Mills equations; the corresponding particle multiplicities are 
given by the relation 

d-^P WtT'^ xO,yO^+QoJ y 

X E 4iPyt{p)ATix)AX{y). (2) 

A=l,2 

Yang-Mills equations with boost invariant CGC initial conditions were solved 
numerically in Ref. [5] to compute the gauge fields at late times and the spectrum 
in Eq. [2] determined. The spectrum matches the pQCD tree level expectation at large 
k± {k± > Qs). Remarkably, for small transverse momenta, the spectrum is best fit by a 
Bose-Einstein form. It has been suggested that such behavior is generic in classical glassy 
non-equilibrium systems |6]. Further, for any finite time, the spectrum is infrared finite. 
One can extract the initial energy density in a heavy ion collision from these numerical 
results, and one obtains e = 0.26 Ag/{g'^As t). Here A^ is a scale of order Q^-estimates 
give As ~ 1.6-2 GeV for RHIC energies. For proper times r = 0.3 fm and g = 2, this 
corresponds to e ~ 20 — 40 GeV/fm^. In addition to providing ab initio estimates for the 
initial energy density and momentum distribution, the classical field simulations provide 
an estimate of the initial eccentricity derived from the universal properties of the glue 
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fields [3- Use of tlie "universal" saturation scale as opposed to other estimates where 
the scale used is process dependent [8] have ramifications for the elliptic flow produced 
in heavy ion collisions; the elliptic flow generated in the Glasma will place a bound on 
the viscosity at RHIC. 

It is important to go beyond these leading order estimates to formulate a 
quantitative theory of the Glasma. A next-to-leading order (NLO) computation 
is important to understand the renormalization and factorization issues (which are 
fundamental to any quantum field theory) and (as we shall discuss shortly) it 
is important to understand the quantum fluctuations which generate the plasma 
instabilities that may speed up thermalization. At NLO in our framework (order 0((7°) 
and all orders in gp), one has two contributions to the average multiplicity [1]. One 
corresponds to the contribution from the production of a pair of gluons while the other 
corresponds to the interference between the LO classical field (in the amplitude) and 
the one loop correction to the classical field (in the complex conjugate amplitude). 
Remarkably, both terms can be computed by solving equations of motion for small 
fluctuations about the classical background field with retarded boundary conditions. 
This means that the average multiplicity at NLO can be formulated as an initial value 
problem and a practical algorithm constructed to compute it. A similar algorithm has 
been constructed and implemented to study quark pair production in the classical CGC 
background field [9]. The gluon pair production NLO contribution is analogous. 



3. The origin of instabilities in the Glasma 

In the Glasma, the classical LO boost invariant fields are purely longitudinal. The 
corresponding momentum distributions (which only generate transverse pressure) are 
very unstable- indeed, they lead to an instability analogous to the Weibel instability 
in electromagnetic plasmas. For a review and relevant references, see Ref. [10]. The 
specific mechanism for a an expanding Glasma is that small rapidity dependent quantum 
fluctuations (that break the boost invariance of the underlying classical fields) grow 
exponentially and generate longitudinal pressure. Romatschke and I performed 3+1- 
dimensional numerical simulations of Yang-Mills equations for a Glasma exploding into 
the vacuum [2]. We chose an initial "seed" spectrum of small rapidity dependent 
fluctuations that are "white noise" Gaussian random fluctuations constructed to satisfy 
Gauss' law. We observed clearly that the maximally unstable modes of the longitudinal 
pressure grow as exp {C\fK^^ with C ~ 0.425-this form of the growth was previously 
predicted for Weibel instabilities in expanding plasmas. The distribution of unstable 
modes also has the "parabolic" shape characteristic of these modes in Hard Thermal 
Loop (HTL) studies of the Weibel instability in finite T QCD. Albeit the 3+1-D solutions 
of the Yang-Mills equations display the same features as those of the HTL studies, a 
deeper understanding of this connection is elusive. 

The trend towards isotropization is suggested by the following observation. At 
early times, the unstable longitudinal momentum mode with the largest frequency 
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(with amplitudes above the background mode distribution) grows hnearly with proper 
time; this corresponds to times when the transverse magnetic fields generated are still 
growing. At the times when these begin to saturate, the maximum unstable mode 
shoots up rapidly. This is accompanied by a decrease in the transverse pressure and 
an increase in the longitudinal pressure. A plausible explanation therefore is that when 
the transverse magnetic fields are very large, the Lorentz force is sufficient to bend 
transverse momentum modes into longitudinal modes. This behavior is seen for initial 
seeds with widely varying amplitudes. The energy density shows a clear deviation from 
free streaming, e ~ r^i'^+^) where 5 ~ 1/16. This value is well below the 5 = 1/3 
required in boost invariant hydrodynamics but is proof in principle that non-trivial 
dynamics can result from rapidity dependent quantum fiuctuations. 

In an interesting recent paper, Fukushima, Gelis and McLerran [1^ computed 
the first quantum corrections to the classical background field at r = 0. The initial 
conditions for the subsequent classical field evolution are quite different from those in 
Ref. [2]. Simulations are underway to determine whether these initial conditions speed 
up thermalization. A full treatment of quantum fluctuations requires that we understand 
the NLO contributions to inclusive gluon production that we discussed previously-these 
go beyond the "WKB" treatment in Ref. [11]. Some of the NLO quantum corrections 
can be absorbed in the evolution of the nuclear wavefunctions with energy while the rest 
contribute to gluon production-a proof of this high energy "factorization" is desirable. A 
related question of great interest is the kinetic theory of the Glasma. How do the Glasma 
flelds decay into particle modes and what are the corresponding kinetic equations ? Work 
in these directions is in progress. 

This work was supported by DOE Contract No. DE-AC02-98CH10886. 
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